SYMMETRY OF BOUND AND ANTIBOUND STATES 
IN THE SEMICLASSICAL LIMIT FOR A GENERAL 
CLASS OF POTENTIALS 

SEMYON DYATLOV AND SUBHROSHEKHAR GHOSH 

Abstract. We consider the Schrodinger operator — h 2 d 2 + V{x) 
on a half-line, where V is a compactly supported potential which 
is positive near the endpoint of its support. We prove that the 
eigenvalues and the purely imaginary resonances are symmetric 
up to an error Ce~ s / h . 



1. Introduction 

In this paper, we study spectral properties of the Schrodinger oper- 
ator 

P(h) = -h 2 3 2 x + V(x) 

defined for x in the half-line (—00, B\. Here h > is the semiclassical 
parameter and V(x) is a piecewise continuous real- valued potential 
supported in [0,5]. 

The operator P(h) with the Neumann boundary condition at B is 
self-adjoint on L 2 (—oo,B); therefore, its resolvent 

R V (X) = (P(h) - A 2 )" 1 , ImA > 0, 

is a bounded operator from L 2 to H 2 for A 2 not in the spectrum of 
P(h). This resolvent can be extended meromorphically as an operator 
L 2 omp — > H 2 oc to A G C with isolated poles of finite rank; these poles 
are called resonances. (The reader is referred to [12] for details.) To 
each resonance A corresponds a resonant state; that is, a nonzero 
u G H 2 oc (—oo, B) solving the equation (P(h) — \ 2 )u = with the Neu- 
mann boundary condition at the right endpoint and with the following 
outgoing condition at —00: 

u(x) = Ae~' lXx l h for all x < and some constant A. 

(Note that for x < 0, u solves the free equation {—h 2 d 2 — \ 2 )u = 0, so 
it must be a linear combination of e ±tXx ^ h .) 
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For Im A > 0, the outgoing condition implies that u is exponentially 
decreasing on the negative half-line and thus «6L 2 ; therefore, A is a 
pole (of the resolvent) lying in the upper half-plane if and only if A 2 is 
an eigenvalue of P(h) on L? . Since P(h) is self-adjoint, all poles in the 
upper half-plane have to lie on the imaginary axis. There may be poles 
A with ImA < and Re A ^ 0; however, we will restrict our attention 
to purely imaginary resonances: 

Definition 1. A positive number k is called a bound state if ik is a 
pole of the resolvent Ry(X), and an antibound state if —ik is a pole. 

We see from above that k is an (anti)bound state if and only if there 
exists a nonzero solution u of the problem 

(P(h) + k 2 )u = 0, (1) 

U X \ X =B = 0, (2) 

hu x ± ku\ x=0 = 0. (3) 

The plus sign in (j3J) corresponds to an antibound state and the minus 
sign corresponds to a bound state. We will also study Neumann eigen- 
values of P(h) on [0,5]; i.e., those k for which there exists a nonzero 
solution to ([T]) with boundary conditions (j2j) and 

u x \ x =o = 0. (4) 

Since the space of solutions to (PQ) and (jSj) is always one dimen- 
sional, bound states, antibound states, and Neumann eigen- 
values never coincide. However, Bindel and Zworski proved in [3] 
that bound and antibound states located away from zero coincide mod- 
ulo errors of order e~ s ^ h for some 5 > 0, if the potential satisfies the 
following conditions: 

3A > 0, V > : V(x) = V for all x e [0, A), 

3e>0: V{x) = for all x G (A, A + e). 

In this paper, we prove a similar result with more general assumptions 
on the potential: 

Theorem 1. Suppose that V is a piecewise continuous real-valued po- 
tential supported in [0, B] and satisfying the following bump condi- 
tion: 

3A > : V(x) > for all x E (0, A). (5) 

Fix two constants < Ck < Ck < oo. Then there exist constants 
C, 5 > 0, ho > such that for h < h and any k e [ck, Ck}: 

1. If k is a Neumann eigenvalue, then there exist a bound state k + 
and an antibound state k_ such that \k — k±\ < Ce~ s ^ h . 



SYMMETRY OF BOUND AND ANTIBOUND STATES 



o 



0.9 



0.8 



0.7 



0.6 



0.5 - 



0.4 - 



0.3 - 



0.2 



0.1 



9©- 



□ B Q u O 



□ 

□ □ 



o 
o 



o 
o 



DD o8 



□ Oo 
□ _ 

□ o □ 

□ ° D n 

□ o □ 

□ no 

□ □ 

o o 
o 



2 3 4 5 6 7 8 



0.9 



0.8 



0.7 



0.6 



0.5 - 



0.4 - 



0.3 - 



o 



> OQ ei 

□ W 

o 



o 



© 



© 



9 



© © 



© 



© 



©- 



o 



2 3 4 5 6 7 8 




-1.5 -1 -0.5 0.5 1 1.5 



-1.5 -1 -0.5 0.5 1 1.5 



Figure 1. Bound and antibound states for two spline 
potentials (splinepot( [0, -0.4, -1, -0.2, -1, 
-0.4, 0], [-2, -1.5, -1, 0, 1, 1.5, 2]) and 
splinepot( [0, 0.2, -1, -0.2, -1, 0.2, 0], [-2, 
1.5, -1, 0, 1, 1.5, 2])) 



2. If k is a bound or an antibound state, then there exists a Neumann 
eigenvalue ko such that \k — ko\ < Ce~ 5 / h . 

The bump condition (JSJ) cannot be disposed of completely, as illus- 
trated by the numerical experiments performed using [2]. Figure 1 
shows two potentials on the whole line, each supported in [—2, 2], and 
the corresponding bound states (denoted by squares) and antibound 
states (denoted by circles). The vertical coordinate of each (anti)bound 
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state on the picture corresponds to its value k; the horizontal coordi- 
nate corresponds to the value of hr x used. We see that the conclusion 
of the theorem does not appear to hold for the potential on the left, 
which does not satify the bump condition; at the same time, it is true 
for the potential on the right. Theorem 1, formulated for the half-line 
case, applies to these numerical experiments on the whole line since for 
even potentials, the set of their (anti)bound states is composed of these 
states for the positive half-line with Dirichlet condition and the same 
states for the Neumann condition; the theorem above can be applied 
with Dirichlet condition in place of ([2]). (However, condition (j3j) cannot 
be replaced by Dirichlet condition in the theorem.) 

The study of resonances in one dimension has a long tradition going 
back to origins of quantum mechanics, see for instance [8]. One of 
the first studies of their distribution was conducted by Regge [TO] ; 
since then, there have been many mathematical results on the topic, 
including [TJ, [I], [5], [6], [TJ, [9], [TTJ, and [13]. Concerning antibound 
states, Hitrik has shown in [6\ that for a positive compactly supported 
potential, there are no antibound states in the semiclassical limit. This 
agrees with our result since there are no bound states in this case. 
Simon proved in [TT] that between any two bound states, there must 
be an odd number of antibound states; the following corollary of this 
result follows almost immediately using the methods we develop to 
prove Theorem 1: 

Theorem 2. Consider the half-line problem with a bounded compactly 
supported potential V (which does not need to satisfy any positivity 
condition). Then for each two bound states < ki < k 2 , the interval 
(k±, k 2 ) contains at least one antibound state. In particular, if there are 
n bound states in some subinterval of (0, oo), then there are at least 
n — 1 antibound states in the same subinterval. 

The proof of Theorem 1 works as follows: we study the evolution 
(in x) of the vectors (u, hu x ) for the three solutions of ([TJ with initial 
data at x = satisfying the conditions ([3]) and (j4]). The idea is to look 
at these three vectors at x — A. Since V(x) + k 2 > on the interval 
(0,A), the transition matrix for the considered vectors from x = to 
x = A will have an expanding and a contracting direction. (In fact, 
if we introduce rescaling x = x/h, then the behavior of the original 
system for small h is similar to the behavior of the rescaled system for 
large x, and the latter will be similar to the behavior of the geodesic 
flow on a two-dimensional manifold of negative curvature.) It turns 
out that our three vectors lie in a certain angle between the expanding 
and the contracting directions, from which it follows that they will 
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stay in this angle for later times (Lemma [2]); what is more, their polar 
angles will get exponentially close to each other (Lemma [7]). Finally, 
we can study how the polar angles of the considered vectors change 
with k (Lemma H]): it follows (Lemma [8]) that the polar angle for the 
solution with Neumann initial data at x = will strictly increase in k 
and the polar angle for the solution with the same data at x = B will 
decrease in k. The proof is then completed by a pertrubation argument 
(Lemma EJ). 

The detailed proofs of Theorems [T] and [2] are given in Section 3. 
Both are elementary and use certain properties of ordinary differential 
equations presented in Section 2. 

The authors would like to thank Maciej Zworski for suggesting the 
problem and for many illuminating discussions. 

2. Preliminaries 

Throughout this section, I is an interval in M., V(x) G L°°(I;M), 
u(x),v(x) G H 2 (I] R), h > 0, and P(h) = -h 2 d 2 x + V(x). Any solution 
to the equation P{h)u = is determined by the vector (w, hu x ) at any 
x; we will sometimes view this vector in polar coordinates: 

Definition 2. Define the length L(u) and the polar angle 8(u) by 
the equations 

u = L{u) cos9(u), 
hu x = L(u) sin8(u). 
Here 9{u) lies in the circle S 1 = M/27rZ. 
Lemma 1. Define the Wronskian W(u,v) by 

W(u, v) = h{uv x — vu x ). 

Then 

W(u, v) = L(u)L(v) sm(6(v) - 6(u)), (6) 
hd x W(u,v) =v ■ P(h)u-u- P(h)v. (7) 

Note that the W(u,v) is just the oriented area of the parallelogram 
spanned by the vectors (u, hu x ) and (v, hv x ). The next lemma tells us 
that if the vector (u, hu x ) falls inside a certain angle in the plane at 
the initial time, then it will stay inside that angle for all later times: 

Lemma 2. Suppose that a 2 < V(x) < b 2 for all x & I and some 
constants a, b > 0. Let u be a solution to P{h)u = and define 

W+(u) = W(u,e bx/h ), W_{u) = W(e- ax/h ,u). 
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Let xq be a point in I and assume that W+(u), W-(u) > at xq. Then 
for x > xq, the functions W±(u) are increasing in x and 

. Uu) 
u > 



VTT¥ 

Proof. We have 



e -bx/h w ( u \ = bu _ hu e ax/h W(u) = au + hu x . 



Therefore, W+(u),W-(u) > yields \hu x \ < bu and thus flSj). Next, 

P(h)e bx/h = e bx/h {V{x) - b 2 ) < 0, 

P(h)e~ ax/h = e- ax/h (V(x) - a 2 ) > 0. 

Using (j7|), we see that d x W± > as long as u > 0. It remains to prove 
that u(x) > for x > x . Suppose this is false and let x\ = mi{x > 
Xq I u(x) < 0}. Then u is not identically zero; since it solves a second 
order linear ODE, L(u) > everywhere. But u > on [xo,xi], so W± 
are increasing on this interval. In particular, W± > at x\ and thus 
(|Hj) holds at this point. However, by the choice of x\ we have u{x\) = 0, 
which contradicts L(u) > 0. □ 

In the next section, we will use the following crude estimate on how 
fast the solutions of an ODE can grow: 

Lemma 3. Assume that |V(x)| < M forx G / and that u is a solution 
to P{h)u = 0. Take Xq,Xi G I; then 

L(u)\ x=Xl < e (i+m*o-*i\/W) . L (u)\ x=X0 . 

Proof. Without loss of generality we may assume that x\ > xq. We 
have L{u) 2 = u 2 + (hu x ) 2 ; thus 

hd x (L(u) 2 ) = 2huu x (l + V(x)) < (1 + M)L{u) 2 

and the lemma is proven by Gronwall's inequality. □ 

Lemma 4. Assume that u(x, k) is a family of solutions to (P(h) + 
k 2 )u = 0, Xq,Xi G /, and u(xo,k) and u x (xo,k) are independent of k. 
Let Q^k) = 8(u(x,k))\ x=Xl , Li(k) = L(u(x, k))\ x=Xl . Then 



2k r xi 



Proof We have W(u,Uk)\ x=Xl = L 1 (fc) 2 0' 1 . (To see that, differentiate 
the formulas in Definition [2] in k and use the definition of the Wron- 
skian.) Now, we differentiate the equation (P(h) + k 2 )u = in k to 
get (P(h) + k 2 )uk = —2ku. It remains to apply (J7J) together with 
W(u,u k )\ x=xo = 0. □ 
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Lemma 5. Assume that $ is a C 1 map from the interval I = [Kq, K±] 
to the circle S 1 = IR/27rZ and $'(fc) > 5 > for all k G I . Suppose 
that ^ : I ^ S 1 is a continuous map such that \^f(k)\ < e < 7T for all 
k. Put v = e/S and I v = [K + u, K\ — u]. Then: 

1. If ko G I v has $(&o) = 0, then there exists k\ G / with = 
^(fci) and \ko — fci| < za 

2. If k\ G has $(&i) = ^/(^i), i/ien inere exists ko G / 
$(/co) = and |fco — fci| < za 

Proof. We can lift $ and \^ to continuous maps : / — > JR.; then 

|^| < e and $(£/) - > <5(Jfe' - k) for jfe' > k. 

1. We have $(&o) = 27rm for some m G Z. Then $(fco + ^) > 
27rm + oV > 27rm + *(& + v) and $(A;o - v) < 2irm + ^{k - u); it 
remains to apply the intermediate value theorem. 

2. Similar to the previous statement, we have = 2ixm + ^>{ki) 
for some m G Z and + z/) > 27rm > — ^). □ 

Lemma 6. Assume that $ is a C 1 map from some interval I to the 
circle S 1 = M/(2ttZ) with > /or a// fc G J. Let * : I -> S 1 6e a 

continuous map such that ty{k) ^ /or a// k <E I. If hi < k% are two 
roots of the equation $ = 0, inen ine interval (/ci,/c 2 ) contains at least 
one root of the equation $ = 

Proof. As in the previous lemma, lift $ and ^ to maps $, \I/ : J — > R; 
we can make < ^(/c) < 27r for all k G /. Since $' > everywhere, 
we have $(%) = 27crrij, where mi < m 2 are some integers. Therefore, 
< 27rm! + ^>(ki) and $(^2) > 27im 1 + ^(^2); it remains to apply 
the intermediate value theorem. □ 

3. Proofs of the theorems 

We assume in this section that < c' k < k < C' k for some constants 
c' k < Ck and C' k > the constants in our estimates will depend on c' k 
and C' k . (We need to expand the interval [c^, C fc ] a little bit to be able 
to apply Lemma O) 

Consider the solutions u±,u ,Ui(x, k) to the equation flTJ in [0, B] 
with the initial data 

w±o(0, k) = 1, d x u (0, k) = 0, hd x u±(0, k) = ±k, 
u t {B,k) = 1, d x u x (B,k) = 0. 

Define Qo(k), Q±(k), and Bi(k) to be the polar angles of vectors 
(u, hu x ) at x = A for u = u , u±, U\. Then k > is 
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• a Neumann eigenvalue if uq and u\ are linearly dependent; that 
is, (recalling that they solve the same second order ODE) if 
2(e (A;) - &i{k)) = 0; 

• a bound state if 2(e+(fc) - 9i(fc)) = 0; 

• an antibound state if 2(Q_(k) — Bi(/c)) = 0. 

Here we count angles modulo 2n. 

To prove Theorem (TJ it suffices to use Lemma [5] (for <£> = 2(9 — ©i) 
and ^ = 2(O — 0±)) together with the following two facts: 

Lemma 7. For some constants C\ and 8\ > independent of h and k, 

\2(e (k) - ± (A;))| < C ie -^ h for all k e [c' k ,C' k \. 

Lemma 8. We have &' (k) - &\(k) > 1/C 2 > for all k G [d k ,C' k ] 
and some constant C% independent of h and k. 



We first prove Lemma [3 Put b = max[ 0] A] V( x ), k b = y/k 2 + b, 
tpo(x) = e~ kx / h , ip + (x) = e kbX ^ h , and consider the Wronskians 

W+(u) = W(u,1>+), W (u) = W(^ ,u). 

These are nonnegative for u = Uo,u± at x — 0. Then by Lemma [21 all 
these six functions are nonnegative and increasing in x for < x < A. 

Our first goal is to get an exponential lower bound on the length 
L(u) for u = Mo, u± at x = A. For w , note that by 

Tl . W^ ,u ) ^ W (uo)\ 

L{uo) ~ - l M * c e ■ 

Same applies to u + . However, w_ needs more careful analysis since 
Wq{uJ) = at x = 0. For that, take < t < 1 and put a = 
mmit A ,A] V{x) > 0, k a = y/k 2 + a, vp^(x) = e ~ kaX/h , and W-(u) = 
W(ip-,u). First, we have by Lemma [3] 

L(u-) > e -a+* a +6)*/(2h) . L(u-)\ X=Q . 
Next, W (w_) > and W+(uJ) > 0, so by © 

W-(u-) > (k a - k)u_tp^ > —L(u-.)if>-. 
Finally, we apply Lemma [2] on the interval [tA, A] to get 

L(u)\ x=A > W -( u ~)\*= tA > L e (Ki-t)-(i+kHb)t)A/h_ 
L(i/)-)\ x= A C 

For t small enough and all k, k(l — t) — (1 + k 2 + b)t > 0, so we have 

i(«-)U > ^ > o. 
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The next step is to use that uq and u± solve the same equation (JTJ) 
and thus W(u ,u±) is constant in x. Therefore, at x = A we have 
by © 

I sm(e(u ± ) - e(u ))\ = ^ Ce ~ kA/h - 

L(u )L(u±) 
That finishes the proof of Lemma [7J 

To prove Lemma [HI first note that by Lemma HI ®[(k) < and 

i r A 

<3>'o(k)> / \u (x,k)\ 2 dx 

ChL(u )%=A Jo 

By (ED, M > L(u )/C. Also, by LemmaEl L(u ) > e c ^- A ^ h L{u )\ x=A 
for < x < A; thus 

£\u (x,k)\ 2 dx > L j\ c ^' h {L{u Q )\ =A ) dx > ^L(u ) 2 \ x=A 

and Lemma [8] is proven, which finishes the proof of Theorem [TJ 

To prove Theorem [2], let $±(fc) = 0(u±)\ x=B ; a bound state corre- 
sponds to 2$ + = and an antibound state corresponds to 2$_ = 0. 
Since 0(u + ) \ x= o is increasing with k, by an argument similar to the proof 
of Lemmalwe get > for all k. Moreover, 2($ + (A;) -$_(£;)) is 

never zero, as this would correspond to u + and w_ being linearly depen- 
dent. We may now apply Lemma[6]with $ = 2$ + and \I/ = 2($ + — $_). 



References 

[1] A. A. Abramov, A. Aslanyan, and E. B. Davies, Bounds on complex eigen- 
values and resonances, J. Phys. A. 34(2001), 57-72. 

[2] D. Bindel, MatScat: MATLAB Codes for ID Potential Scattering, 
http : / / cims .nyu . edu/~dbindel/matscat 

[3] D. Bindel and M. Zworski, Symmetry of Bound and Antibound States in the 
Semiclassical Limit, Lett. Math. Phys. 30(2007), 10-20. 

[4] P. Briet, J.-M. Combes, and P. Duclos, On the Location of Resonances in the 
Semi- classical Limit LL, Comm. Math. Phys. 12(1987), 201-222. 

[5] R. Froese, Asymptotic distribution of resonances in one dimension, J. of Dif- 
ferential Equations 137(2), 1997, 251-272. 

[6] M. Hitrik, Bounds on scattering poles in one dimension, Comm. Math. Phys, 
208(1999), 381-411. 

[7] E. Korotyaev, Inverse resonance scattering on the real line, Inverse Problems 
21(2005), 325-341. 

[8] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic The- 
ory, Third Edition, Elsevier, 1977. 

[9] L. Nedelec, Asymptotics of resonances for matrix valued Schrodinger opera- 
tors, |math.SP/0509391[ 
[10] T. Regge, Analytic Properties of the Scattering Matrix, Nuovo Cimento 
10(1958), 671-679. 



10 



S. DYATLOV AND S. GHOSH 



[11] B. Simon, Resonances in One Dimension and Fredholm Determinants, J. 

Fund. Anal. 178(2000), 396-420. 
[12] S.-H. Tang and M. Zworski, Potential Scattering on the Real Line, lecture 

notes, http : //math.berkeley . edu/~zworski/tzl . pdf 
[13] M. Zworski, Distribution of poles for scattering on the real line, J. Funct. 

Anal. 73(1987), 277-296. 
E-mail address: dyatlov@matri.berkeley.edu 

E-mail address: subb.ro@matli.berkeley.edu 

Mathematics Department, University of California, Evans Hall, 
Berkeley, CA 94720, USA 



